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ON THE NODAL LINES OF EISENSTEIN SERIES ON 
SCHOTTKY SURFACES 

DMITRY JAKOBSON AND FREDERIC NAUD 


Abstract. On convex co-compact hyperbolic surfaces X = r\H 2 , 
we investigate the behavior of nodal curves of real valued Eisenstein 
series F\(z, £), where A is the spectral parameter, £ the direction at 
infinity. Eisenstein series are (non-L 2 ) eigenfunctions of the Lapla- 
cian Ax satisfying A xF\ = (| + X 2 )F\. As A goes to infinity (the 
high energy limit), we show that, for generic £, the number of in¬ 
tersections of nodal lines with any compact segment of geodesic 
grows like A, up to multiplicative constants. Applications to the 
number of nodal domains inside the convex core of the surface are 
then derived. 


1. Introduction 

Let HI 2 be the hyperbolic plane endowed with the usual metric of con¬ 
stant negative curvature —1. Assume that T is a convex co-compact 
group of isometries, i.e. a Schottky group with no parabolic elements, 
and denote by X = r\HI 2 the quotient surface. Such a surface has 
infinite area and the ends are hyperbolic funnels. Let Ax denote the 
hyperbolic Laplacian on A". Its L 2 -spectrum has been described com¬ 
pletely by Lax and Phillips in [12]. The half line [1/4, -|-oo) is the 
continuous spectrum and it contains no embedded eigenvalues. Let 
<5(r) be the Hausdorff dimension of the limit set A(r) of T. The limit 
set A(T) is defined has the set of accumulation points in dM 2 of the 
orbit of any point z G HI 2 under the action of T: 

A(T) := T\z fl <9HI 2 . 

The rest of the spectrum (point spectrum) is empty if S < |, finite 
and starting at 5(1 — 5) if 5 > The fact that the bottom of the 
spectrum is related to the dimension 5 is due to Patterson [15]. One 
way to parametrize the continuous spectrum is through the so-called 
Eisenstein Series. Before we can give a formal definition of Eisenstein 
series, let us recall that under the above assumptions, when T is non- 
elementary (i.e. X is not a hyperbolic cylinder), A" can be decomposed 
as 

A = Xq U 3q U ... U T n/ , 
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where X 0 is a compact surface with geodesic boundary and 3b,..., 3 r n/ 
are the funnels. Each funnel 3y is isometric to the cylinder 

(0, 2]p x (WL/£jZ)g, 

endowed with the conformally compact metric 

A 2 _ dp 2 + (1 + p 2 /4)dd 2 

where p = 0 corresponds to infinity and £j is the length of the geodesic 
boundary at p = 2. Let Rx(s\z,w ) denote the Schwarz kernel of 
the resolvent (Ax — s(l — s)) _1 which by Mazzeo-Melrose [BE] has a 
meromorphic continuation (in s ) to the whole complex plane. Then if 
s is not a pole, the limit (using the above coordinates in the funnel) 

Es(z,Q ■■= lim p~ s Rx(s;z,(p,0) 

p—t o 

exists and dehnes an eigenfunction of the Laplacian A \E s (z,f) = s(l — 
s)E(s; z,£), parametrized by a point £ at infinity, called Eisenstein 
series. In particular if s = 1/2 + *A, we have 

A x E,(z,()= fj + A 2 Ws,f)- 

These Eisenstein series, like their analog in the finite volume case, 
provide an explicit spectral resolution of the Laplacian [2]: 

2\dn x (\,z,z') = ICW^ + AH 2 y. f 1 ’ E 1/2+a (z,Z)E 1/2 _ a (z',Z)<K, 

Z7T j.1 J o 

where 

c(s] = 2 ~‘ r < s ) 

v/ir(s- 1/2)' 

In this paper we want to investigate the zeros sets of high energy Eisen¬ 
stein series , so we consider real valued Eisenstein functions i.e. we take 
real parts (which are again eigenfunctions) and set for all z G X and £ 
a direction at infinity, 


AM) :=Re (Si +a M)). 

Below we show a plot in the Poincare disc model for a symmetric two 
generator Schottky group with £ = i and A = 30. 
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It is a natural question to investigate the shape and behaviour of the 
zeros sets (also called nodal lines in dimension 2) of F x (z ,£) as the 
frequency A goes to infinity. For genuine L 2 -eigenfunctions on compact 
manifolds, there is a tremendous amount of work in that direction, 
and we refer the reader to the recent survey [18]. However, in the non¬ 
compact case and infinite volume case, this seems, to our knowledge, to 
be the very first related work. Numerical experiments show that nodal 
lines exhibit a mixed behaviour: horocyclic shape close to infinity (as 
depicted in the above picture) while in the compact core they look 
more like a genuine high energy eigenfunction, we refer the reader to 
§5 for a high energy plot with A = 150. 


Even in the case when T is elementary, the numerics show a highly 
non trivial nodal structure. Below we plot the Eisenstein series F x {z, £), 
in the Poincare half-plane for £ = 2.5, A = 40. The group is generated 
by z i —y e l z with £ = 1.5. 
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In the case when <5(T) < 1/2, then the lift to HI 2 of F\(z,£) admits 
a convergent series expression (in the unit disc model), see [5], Lemma 
5, for a proof of that fact. 



1/2 

COS 





) 


where z6i 2 and £ G <9HI 2 belongs to the domain of discontinuity of Y 
that is <9HI 2 \ A(r). For each term in this sum, the phase function has 
its level sets on Horocycles based at the point 7 -1 £ at infinity, so it is 
basically a superposition of hyperbolic plane waves. 

Because F\ is an eigenfunction of an elliptic operator with real an¬ 
alytic coefficients (the hyperbolic laplacian), it is automatically a real 
analytic function. The nodal sets Na(£) are defined as usual by 


X a (£):={zGX : F A (z,£) = 0}. 


These sets are real analytic curves (with possible isolated singular 
points) and therefore rectifiable. Let a denote the length measure 
induced on N A (£), then by translating almost verbatim the arguments 
of Donnelly-Feffermann [5] (which is a purely local proof ), one ob¬ 
tains that for all compact K C X with non-empty interior, there exists 
Ck > 0 such that as A —» oo, 

C^X<a(X x (OnK)<C K X. 


In this paper we go beyond by proving the following result. Given a 
geodesic C, we will define a notion of £-non symmetry (£-NS), see §3, 
which rules out cases where the geodesic 6 is an axis of symmetry for 
certain geodesics related to £. The following holds. 


Theorem 1.1. Assume that <5(T) < Let C be a geodesic which 
satisfies £ '-NS. Then for all compact non empty segment Co C G, one 
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can find a constant Cq such that as A goes to infinity, we have 

Co 1 a < #(k a (0 n e 0 ) < c 0 a. 

The above statement is non-empty : for all geodesic C, £-NS is sat¬ 
isfied for almost all directions £ at infinity, see §3. The upper bound is 
actually valid in greater generality for real analytic curves and generic 
£, see comments in §3. 

We point out that several recent papers also focus on proving upper 
and lower bounds on the number of intersections of nodal lines with 
geodesics segments. On compact non positively curved surfaces with 
boundary, Jung and Zelditch m , show that #(N A nC 0 ) goes to infinity 
as A goes to infinity, when C is a boundary curve. On the other hand, 
a similar statement holds [TQ] on a negatively curved surface (without 
boundary) and when G satisfies a non symmetry condition. On the 
modular surface PSL 2 (Z)\IHI 2 , Jung [9] obtains effective lower bounds 
of the type 

c^Af -6 < #(\ne 0 ), 

for the Maass-Hecke eigenfunctions (with discrete spectral parameter 
A*, as in our case) for a large portion of A/Js and when 6 is a vertical 
geodesic segment in the modular domain. In [7], Ghosh, Reznikov and 
Sarnak, assuming Lindelof’s hypothesis, obtain a related lower bound 

^- £ <#(x A ne 0 ), 

for all Afc large enough. On the Flat 2-torus, Bourgain and Rudnick [3] 
were able to show that for non geodesic curves, 

#(^a n e 0 ) > ca 1-£ . 


It seems to us that Theorem 1.1 is the only optimal counting result 
so far. Of course our setup of infinite volume is helping us somehow, al¬ 
though there are some different technical difficulties to overcome. The¬ 
orem [Li] is a consequence of the following (restriction) equidistribution 
result, which is of interest in itself. 


Theorem 1.2. Let T be a convex co-compact group with J(T) < |. 
Let Go be a finite length geodesic segment of a geodesic satisfying £ -NS. 
Then for all <p G G 1 (C 0 ), we have 

lim [ (F x (x,£)) 2 ip(x)da(x) = \ [ E 1 (x,f)<p(x)dcr(x), 

A—H-oc J Cq J Co 

where is the positive harmonic Eisenstein series at s = 1, and 

a stands for the length measure. More generally, the same statement 
holds for all real-analytic compact curve C 0 , for a generic choice of £. 


This above theorem is a ’’restriction” version of the main equidistri¬ 
bution result of [8j, and this is where the £-non symmetry assumption 
is required. Similar equidistribution restriction results are known on 
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compact manifolds (so-called ”QER” ) when the geodesic flow is er- 
godic, also under a non symmetry assumption, see for example Toth- 
Zclditch HZ). We also refer to the paper of Dyatlov-Zworski [6] for 
a semi-classical framework that generalizes the preceding results. See 
also Bourgain-Rudnick [3J for related results on the torus. 


Most of the above mentioned works are motivated by the study of 
nodal domains. It is a notoriously challenging problem to count them 
and dHDUE] provide the very Erst (deterministic) examples of eigen¬ 
functions where one is actually able to show that the number of nodal 
domains goes to infinity at high frequency. As a corollary of theorem 


A"o denote the convex core of X (the compact part with funnels re¬ 
moved) and let M%( A) be the number of (open) connected components 
of 


1.1 we prove the following. Assume that T is non-clcmentary, and let 


Int(Xo) \ tt A (0. 


Corollary 1.3. Under the above hypotheses, for almost all f, there 
exists a constant C > 0 such that for all A > 1, we have 

Aff(A) < CX 2 . 

If T is elementary i.e. r\HI 2 is a hyperbolic cylinder, let Co denote 
the unique closed geodesic in X = r\HI 2 , and let C(r) be the collar of 
size r > 0: 

C(r) := {z G X : dist(z, C 0 ) < r}, 

and let M^(A) denote again the number of (open) connected compo¬ 
nents of 

Int(C(r))\N A (0. 


Corollary 1.4. Using the above notations, for almost all f, there exists 
a constant C > 0 such that for all X > 1, we have 

Mt(X) < CX 2 . 

It is important to notice that these upper bounds, which are analogs 
of Courant’s nodal domain theorem, are not obvious facts: eigenfunc¬ 
tions F\(z,£) do not satisfy any boundary condition on dX 0 or dC (r). 
It is tempting to believe that this bound is optimal, but we have no 
serious clue so far. 

The plan of the paper is as follows. In §2 we recall some basic facts 
about hyperbolic planes waves and Eisenstein Series. In §3 we prove 
Theorem |1.2| and a result on the asymptotic average on a geodesic seg¬ 
ment (Proposition |3.1|). Theorem 1.2 will be used for both lower and 
upper bounds in the proof of Theorem 1.1 while Proposition 3.1 is 
critical for the lower bound, see §4 for details. We point out that while 
the lower bounds and the equidistribution result rely on elementary 
real analysis (stationary and non-stationnary phase principles for os¬ 
cillatory integrals), the upper bound requires some complex analysis. 
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This problem is already present for compact problems where the upper 
bound of Donnelly-Feffermann has not yet been proved in the C°° cat¬ 
egory (Yau’s conjecture). Because the methods we use here are fairly 
elementary and robust, we expect these set of results to be extendable 
to variable curvature cases, with a negative pressure condition, as long 
as some analyticity is available. 

Acknowledgments. This work was mostly done while FN was a 
member of UMI 3457 at universite de Montreal, supported by CNRS 
funding. Both authors are supported by ANR ”blanc” GeRaSic. DJ is 
also supported by NSERC, FQRNT and Peter Redpath Fellowship. 


2. Basic estimates and convergence 

In this section we gather various basic estimates that wil be required 
later on. We start with some facts on Busemann functions that will be 
used frequently throughout the paper. 


2.1. Busemann functions. We will mostly work with the unit disc 
model 

H 2 = D = {z G C : |z| < 1}, 
endowed with the metric 

, 9 4 dzdz 


The hyperbolic distance between 0 and z is given by 
(1) d(0,z) = log Q + . 

Given two points z,w G H 2 and £ G <91HI 2 , i.e. |£| = 1, the Busemann 
function B^(z,w) by 

Bt(z,w) = lim d(z,£ t ) ~ d(w,£ t ), 

t—y+oo 


where t i-A £* is converging to £ as t -A +oo. From that definition ones 
deduces several standard properties of Busemann functions which can 
be checked easily. 

• For all z : w,y , we have B^(z,y) = Bg(z,w) + B^(w,y). 

• For all z,w, B^(z,w) = —B^(w,z). 

• For all isometry g of the hyperbolic plane, we have B^(gz, gw) = 

Bg-l^Z.W). 

• The formula holds: 


Bt(z,w) = log 


(1- 

\w 

l 2 )l^-el 2 

(1- 

\z I 

2 )k-£l 2 
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The level sets of ^ i— B^(z,w) are Horocycles based at £ . The hyper¬ 
bolic analog of monochromatic plane waves are functions of the form: 

* e iXB ^°’ z \ 


It is shown in [8] that if 5(T) < I, then generalized eigenfunctions 
Ei/2+i\(z , £) which are a priori defined through analytic continuation 
admit the convergent series formula 


z 7er 


e (|+*A)B e ( 0 , 7 z) 


In particular, the formula for the (real) Eisenstein series becomes 

Fx(z,0 = ^e^^cos(A^(0, T ^)). 

7 er 

We start by a simple estimate. 


Lemma 2.1. Assume that z E K, where K is a compact subset of 
HI 2 . Then for all multi index a = (an ,... ,chn), there exists a constant 
C ( K , a) such that 

\d a | < e~2 d ^C(K,a), 

where z = x i + ix 2 and d a = ... » 9 . 

OXoc-^ OXotjy 

Proof. We only compute the first derivatives, the rest follows by an 
easy induction. Writing 

_ e |(B e ( 0 , 70 )+B 7 _i e ( 0 ,z)) 

we have for j = 1 , 2 : 

g fi^(o,7 = ( 2x i _ 2 ((7~ 1 0j ~ x j) \ §b £ (o, 7 o) 

J V ) V l-kl 2 | 7 _ 1 £--I 2 ) 

Now remark that by formula (jl| we have 

B,( 0 . 7 0 ) = log (7)1) = —<i( 0 , 7 0 ) + 2 log (7W7) . 

Because £ is not in the limit set of T, the distance I 7 O — £| is uniformly 
bounded from below, so there exists a constant C\ > 0 such that 

B <:{ 0,70) < -d(0, 7 0) + C' 1 . 

Since z is confined to a compact set K and y _1 £ remains on the unit 
circle, the proof is done. □ 

This simple estimate implies that if £(T) < I, then the series defin¬ 
ing F\(z, £) (and the derivatives) are uniformly convergent on every 
compact subset of HI 2 . Indeed we recall that Poincare Series 

P r (s) = ^V sd(0 ’ 70) 
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are convergent for all s > 5(T), see |J5]. Therefore, F\(z ,is a C°° 
function, which is not surprising. From this elementary estimate, we 
have readily the following consequence which is worth highlighting: 
for every compact set K C HI 2 , there exist Ck > 0 and Ck > 0 
independent of A >> 1 such that 

(2) \\F\\\l°°(k) < Ck- 

Notice that in the compact or finite volume case, L°° norms of high 
energy eigenfunctions are usually not expected to be bounded: For 
example Maass wave forms (L 2 eigenfunctions) on the modular surface 
PSL 2 {fL)\ HI 2 are not L°° bounded, see [ IB] . 


3. Restriction Theorems 


In this section we prove the main equidistribution theorem for re¬ 
striction to geodesics (and more) stated in the introduction. Most of 
the results will rest on repeated applications of stationary and non¬ 
stationary phase formulas. All the computations will be done in the 
disc model, but of course results do not depend on the choice of a 
particular model for HI 2 . 


3.1. Asymptotic average on a geodesic segment. Geodesics in 
the disc model will be parametrized in the following way. Let g be a 
Moebius map of the unit disc, then the image of 

g : (—1, +1) —>• HI 2 

is a geodesic. We denote it by G g . Conversely all geodesics of the disc 
can be viewed that way: given a geodesic C and a point z 0 G C, there 
exists a Moebius transform g such that g(0) = z 0 and g((— 1, +1)) = C. 
What we first prove is the following. 

Proposition 3.1. Assume that 5{T) < LetG ro be a geodesic segment 
in HI 2 , parametrized by g : [—r 0 , +r 0 ] — > HI 2 . Then there exists a non 
empty open interval J C [—r 0 , +r 0 ] and C > 0 such that as A goes to 
infinity, we have 

C 

sup / F\(g(r),f)dr 

a</3£j J a 



Proof. Since 5 < |, using the representation of F\ as a sum of conver¬ 
gent series we are left with estimating the sum 



e (|+a)B e (0,7g(r))^, r 


where a < f3 G J C [—To, +r 0 ] and J has to be chosen. The choice of J 
will follow from a careful analysis of the stationary points of the phase 
B^(0, g'f(z)). Writing 

A;(0,7 S ( r )) = £*(0,70(0)) + £5-17-1(0(0, r )> 
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we deduce that 

(3) |W0.7 g ( r))) = 2 (1 _^-^ + ^ , 

where we have set a 7 = Re(r/ H 0 )> b i = Im (# l i HO)- The 
critical points are then given by 

7,7 = ^ i 1 ± a / 1 ~~ ° 7 ) 

if a 7 7^ 0, and 0 otherwise. Remark that only r~ can be a critical point 
(r 7 is outside the disc) and we have if a 7 7^ 0 , 

^( 5 H°> 7 ^C r ))) > ||a 7 ||r-r+||r : ^r“| 

> ||r - r 7 |(1 - r 0 ). 

More precisely, consider the continuous, injective map 

F : [—1, +1] -A [—1, +1] 

defined by 

_. , x 

F ( x ) = , , n - 


v ' i + ^r^’ 

then F(a 7 ) is the unique possible critical point of the phase. In all 
cases, we have a lower bound for the derivative: for all r G [ — r 0 , +r 0 ], 

7(B { (0, 7 j(r))) >C(r„)|r-F(aJ|i 

where C(r 0 ) is uniform in 7 . The goal is now to find a non empty 
interval that is uniformly away from all the critical points. Let us 
consider _ 

K:= U^" 1 °7- 1 (0 COT 2 , 

7 er 

and set for all z G cffl 2 , F(z) = F(Re(z)). Then dehne 

3 := F(K) fl [—r 0 ,+r 0 ], 

then “B is a compact subset of [—r 0 , +ro] which contains all the possible 
critical points of the phases. Observe now that because F is injective 
and continuous there exists rj > 0 such that 

3 = F(K \ 77 ) U _D(+1, 77 )), 

where D(z,r /) := {|w| = 1 : \z — w\ < 7 }. Since F is smooth away 
from —1 and + 1 , we deduce that 

dim^(B) < dim h(K), 

where dim ^ stands for the Hausdorff dimension. Because we have 
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and since the set of accumulation points of the orbit T.£ is exactly the 
limit set A(r), we deduce that 

dim h (3) < S(T) < 1 . 

As a consequence, [—ro, +ro] \ B has non empty interior. We therefore 
pick J C [—r 0 , +r 0 ] an interval such that 

1 ns = 0. 


On this interval J all points are uniformly away from the ’’bad critical 
set” B. This will allow us to use the following version of non-stationary 
phase estimate. 


Lemma 3.2. Let I be a compact interval. Let <f> G C 2 (I) and tp G 
C l (I). Assume that for all x G I,&(x) 0. Then one can find a 

constant M(I ) such that for all a < b G I, for all A > 1 we have 

< max(l, l|$ // ||c°(r))ll^||ci(J) 

A(inf/Id *'!) 2 



The proof of this fact is elementary: just integrate by parts. We now 
apply the above non-stationary principle to each term in the sum 

f e (|+a)s ? (o,7 9(r))^ r 

7er 'a 


We use the fact that for all r G J, 


d 


— (B c (0,7t/(r))) 


> C(J) > 0 , 


and apply Lemma 2.1 to deduce that uniformly in 7 , 


rh 1 

e (2+*A)B c (0,7(/(r))^, r _ q 


g2 d (°D0) 

A 


One has also to check that || r yg(r)))\\ c o^ is uniformly bounded 

from above, which follows easily from the formula (J3]) . The end of the 
proof follows from convergence of Poincare Series. □ 


3.2. The condition of £-non symmetry. We will first state the 
definition on the universal cover. Given two different points 7 ^ r / 2 G 
S 1 := <9HI 2 , we denote by Cthe unique (non oriented) geodesic in 
H 2 whose endpoints are 71 , 1 / 2 - 

Definition 3.3. Lett; G dEI 2 \A(r). LetG g := C 9 (_i) j5 ( + i) = < 7 ([—1, +1]) 
be a parametrized geodesic as above. We say that G g is f -non symmetric 
(f-NS) iff 

(1) V 71 7 ^ 72 G T, G g and C 7l £, 72 £ are non orthogonal. 

(2) V 71 7 ^ 72 G T, G g 7 ^ 67 ^, 72 ^- 
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Remark. Condition (1) implies that for all 71 7 ^ 72 , we have 
Re(c / _1 71 (0 - c/ _ 1 7i(0) ^ 0. 

Indeed, we cannot have g i_ 1 7 i(£) = 9 _ 1 7i( 6 otherwise we would have 

7b 1 0 72(0 = 6 

which is impossible outside the limit set (recall that £ ^ A(T)). There¬ 
fore 


Re (.9 ' 71 ( 6-9 Si( 6 ) = 0 =» 9 _ 1 (7i0) = 9 
which by conformal invariance of angles implies that 

Gg _L 

Note that condition ( 2 ) is always fulfilled if G g is a trapped geodesic , 
i.e. both endpoints g(l) and g(— 1) belong to the limit set A(T). We 
prove below that these conditions have full measure with respect to £. 

Proposition 3.4. Let G be a geodesic. Then for Lebesgue almost all 
£ G S 1 \ A(T), C satisfies f-NS. 

Proof We assume that G is parametrized by 

C — g{[~ 1) +1]); 

for some Moebius map g. First we remark that either g( 1) belongs to 
A(r) and (2) is automatically satisfied or g{ 1) ^ A(r) and its orbit 
under the action of T is discrete in S 1 \ A(r). Therefore, ( 2 ) is satisfied 
if we chose f to belong to 

s 1 \ ^A(r)u y{7j(l)}J , 

which is a set of full measure in S 1 \ A(T). 

If (1) is violated for 71,72 G T, we must have 


9 _ 1 (7i0) = 9 1 (726). 

This identity can hold for only finitely many £ e S 1 . Indeed if h 1 ,h 2 
are two (orientation preserving) isometries of the hyperbolic disc, the 
equation 


(4) MO = * 2(0 

has at most two solutions in S 1 = <91HI 2 : for orientation reasons this 
equality cannot hold identically on S 1 , any solution of Q is a root 
of a non zero polynomial of degree at most 2. We therefore have to 
remove from S 1 \ A(T) a countable set of possible solutions to make 
sure that (1) is satisfied. I 11 a nutshell, both (1) and (2) are satisfied 
for all £ G S 1 \ A(T) except for a countable set, the proof is done. □ 
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On the quotient surface X = T\HI 2 , the condition £-NS translates as 
follows. Given a geodesic C on the surface, it satisfies £-NS if C is never 
equal or orthogonal to geodesics that start and end at £ (at infinity). 
Indeed, geodesics that start and end at £ are lifted on H 2 to geodesics 
whose endpoints are equal to £, mod T, that is geodesics of type C 71 ^ i72 ^, 
for some 7i 72 € T. 


3.3. Proof of the equidistribution result on geodesics. The goal 
of this subsection is to prove the following fact, which implies straight¬ 


forwardly Theorem 1.2 


Theorem 3.5. Assume that V is a convex co-compact group with <5(T) < 
|. Let G = g([— 1,+1]) be a geodesic satisfying (-NS. Then for all 
0 < r 0 < 1, for all <p E C^Q-ro, +r 0 ]), 


»+ro 


'■+ro 


lim 

A—^--(-oo 


(F\(g(r),0) <p{r)dr = \ / Efigir), £)(p(r)dr. 


-ro 


'-r 0 


Proof. We start by writing 

(F x (z,0) 2 = l\E 1/2+iX (z,0\ 2 + ^Re ((E 1/2 +jA (^,0) 2 ) , 

so that we have to investigate 


, >+ro 


"+ro 


(F X {g{r),Q) (p(r)dr = \ / \E 1/2+iX (g(r), £)\ 2 ip(r)dr 


' —r 0 


r -ro 


7(A) 


I Re 


r+r 0 


( E 1/2+iX (g(r),€)) 2 (p(r)dr 


' —r 0 


7(A) 

We will first analyze I i(A). By uniform convergence we can write 

/ + r 0 l 

e 2( S d0,7i9W)+^(0,729(r))) e a$ 71 , 72 (r)^( r ^ r; 

-rn 


where 


7 l, 72 Sr J ~ r ° 

$ 7i,7 2 ( r ) = ^(°,7i 9(r)) - Be(0,j 2 g(r)). 


Writing 

$ 71 , 72(0 = 7W(0))-R0°> 723(0))+R 9 -i 7i -i ? (0, r)-B g - l7 -^(0, r), 

we deduce that 

d (r)) _ 9 R e 0 _1 72 _ 0 - g-^Oir 2 ~ 1) 

1^71,72 V )) ~ Z 


dr 

and therefore, 

inf 

[-r 0 ,+r 0 ] 


\r-g~ 1 7 2 0l 2 l r -wOlOl 2 


^($71,72(0) 


> 


1 - 0 


16 


|R e(g x 7 2 1 £-9 OOOI 
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Because we are assuming property £-NS, part (1), we know that for all 
71 7 ^ 72 this lower bound cannot vanish. This will allow us to apply 
the non-stationnary phase Lemma 3.2 to the off-diagonal sums above. 
More precisely, we have 

r>+ro 


E 


/ -\-ro 

e( B «(0’ 79 ( r) V(r)dr 

r° 

r, ~H *0 1 

e 2( B ((°’'n9{r))+B i (0,'y 2 g(r))) e i\^ in2 (r) 


71^72 70 

where we can write again by uniform convergence 


E 

7 er ■ 


»+ro 


r+ro 


e (s e (o, 7 ' 9 ( r ))(p(r)d r = / Ei(g(r), £)(p(r)dr. 


-ro 


' —ro 


ft is important to notice that Ei(z,£) is a positive non vanishing Har¬ 
monic function on the unit disc which satisfy the trivial lower bound 
(given by the identity term in the sum): 


To complete the asymptotic analysis of /i(A), we therefore have to show 
that the off-diagonal contribution goes to zero as A goes to infinity. Let 
us write 

_ r+ro | _ 

/ e 2^(«4i9W)+B £ (0,72 9 W)) e «A$ 71 , 12 (r)^ r ) (Ir = ^ /^(A). 


71^72 


' —ro 


717^72 


By Lemma 2.1, we have uniformly in A, 


(5) 


7i?72 


(A) | < C'(r 0 )e _ 2 d ( 0 ’ 7 l 0 ) _ 2 d ( 0 ’ 72 °) 


while by the above analysis of phases $ 7l 72 and Lemma 3.2, we do 
have for all 71 7 ^ 72 , 


( 6 ) 

Because we have 

we can deduce that 


IW*)l=o(jV 


^ e -2 d (°.Tl 0 )-2 d ( 0 .72°) < +OG 


71 >72 


iim y; WA) = 0. 

A —^+00 £ 


717^72 


indeed, fix e > 0, and choose T so large that 


C(r„) x 


E 


3 — 2 rf (°771°)-2 d (°’720) 


71^72 

d(0,7i0)>T or d(0,720)>T 


~ 2’ 
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where (7(Vo) is the constant in estimate (J5]) . Writing 


'y i -^ 71,72 (-^) 

71^72 


e 

— 2 


E 


n 7 ii 72 ("v I) 


717^72 

d(0,7^0)<T and d(0,720)<T 


using i§. we can choose Aq so large that for all A with A > Aq, 


E 


1^71:72 (^) I — 2’ 


71^72 

d(0,7i0)<T and d( 0,720)<T 

and we are done. 

Next we move on to the analysis of / 2 (A). Again using uniform 
convergence, we have 
r+ r o 

I (E 1/2+iX {g(r),€)) <p(r)dr= , %2 

■^71.72 (^)> 

71172 


' —r 0 


where 


^71)72 (“^) — 


r~\~ r 0 1 

/ e 2( S d 0 >7l9h))+^(°-72ff(r))) e iAe 71i72 (r)^|- r ^ r ^ 

' -r 0 


with 


® 7i)72 (0 = ^(0,7i 9(r)) + £ € (O,720(r)). 

Using the same tricks as above, one can compute 

d _ r 2 a 71 — 2r + a 71 r 2 a 72 — 2r + a. 


+ 2 - 


'72 


dr (©Tl.TaW) 2 ^_ r 2 )| r _^-i 7i i ^|2 1 “(1-r 2 )|r - g~ x ^ X £| 2 


-^ 71,72 (0 


1 - r 2 |r - cr^ 1 ^| 2 |r - ^| 2 ’ 


where 


P 7i7 /(r) = (a 7 + a 7 /)r 4 — 4(a 7 a 7 / + l)r 3 + 2(a 7 + a 7 /)r 2 
—4(a 7 a 7 / 4- l)r -1- o 7 -1- oy, 

and a 7 = Re(<? _1 7 _1 £). Therefore we get the lower bound 


M) 


> -LR 


71:72 


:oi. 


A key observation is that this polynomial has always degree 3 or 4. 
Indeed, if we have 

<z 7l T a 72 — 0 and o 7 l ci 72 -f 1 = 0 , 
then (a 7l , a 72 ) G {( 1 , — 1 ); (— 1 , 1 )}, which would mean that either 

7 i _1 £ = g(~ 1 ), 7 2 _1 £ = 9( 1 ) 

or 

7i _1 £ = 0 ( 1 ), 7 2 _1 £ = 9(- !)■ 
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This is not possible because of condition (2) in £-NS. To conclude the 
proof, we will need the following Van der Corput’s style Lemma, to 
deal with the possibly highly degenerated stationary phases. 

Lemma 3.6. Let I be a compact non trivial interval and F e C 2 (I), 
tp G C' 1 (J). Assume that for all x E I, we have 

\F'(x)\>C\P{x)\, 

where P(x) is a polynomial of degree d > — oo. Then as X goes to 
infinity, we have 

e lXF[x) p{x)dx = o(x~^). 



Proof Let P(x) = ao+aix+...+adX d , with ad ^ 0. Let x\, X 2 ,... ,Xd G 
C be the roots of P(x) so that we can write 

(7) P(x) = a d (x-x 1 )...(x-x d ). 

Let e > 0 to be specified later on. For all e > 0 small enough, set 

/(e) := {i 6 / : Wj = l,...,d, \x — Xj\>e}. 

Then for all e > 0 small enough I e is a finite union of closed intervals 

d 1 

1(e) = |J /,(£), 

t= 1 


with d' < d independent of e. On each interval /g(e), F' does not vanish 
so that we can integrate by parts 

ip(x) 


e iXF W<p(x)dx = -X- 


'L(e) 


1 

e i\F(x) 

iX 

d 

( <p(x) 

dx 

\F'(x) 

x E /(e), 


JhT) 


\F\x)\>C\a d \e d , 

which yields for all A > 1 and all e small, 


F'{x) 
dx. 


9I t (e) 


J\F(x) 


'm 


(p(x)dx 


< 


C 

A^7’ 


where C is independent of A,e. Writing 


e iXF ^<p{x)dx = I e lXF(x) (p(x)dx + 
Jl(e) 

= 0(e) + O 


i\F(x) 


h\I(e) 


e iXF{x) (p(x)dx 


Xe 2d J ’ 


1 

e = A 2d + 1 , 


we then choose 
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and the proof is done. □ 


Note that the rate of decay as estimated above is far from being 
optimal, bnt enough for our purpose. We can now finish the proof of 


the equidistribution theorem. By Lemma 2.1, we have uniformly in A, 


|J 71 i72 (A)| < C'(r 0 )e“^ (0 ’ 7l0) -^ (0 ’ 720) 


3.6 


while Lemma 
individually as A goes to +oo, 


and the computation of 0' (r) above show that 


|J 71 , 1 i (A)| = o(aA) 


The same arguments as above then yield 

lim J 2 (A) = 0, 

A—^-|-oo 


finishing the proof of Theorem 3.5 □ 


3.4. Equidistribution on real analytic curves. In this section, 
we explain in a nutshell how the above equidistribution theorem on 
geodesics can be extended to all real analytic curves, for almost all £. 
The ideas are very similar to the above proof, but the price to pay to 
obtain a result at this level of generality is that the generic conditions 
on £ have no longer a simple geometric interpretation as in the £-NS 
statement. We have chosen to include details on this generalization 
because it could be useful in some situations. 


Without loss of generality, we will assume that g is a Moebius map 
of the unit disc and that i : [—r 0 , +r 0 ] —* HI 2 is a real analytic complex 
valued map with f(0) = 0 and £'(r) ^ o for all r e [—r 0 , +r 0 ]. We will 
consider the map 

got: [-7-0, +r 0 ] ->■ H 2 

as a parametrized curve on which we want to prove the same statement 
as above. Following the exact same lines, we need to analyze the two 
phase functions 


$7i,72( r ) = £*(°,7is00))) - 72^00))), 

©71,72 0) = 5 0°>7W00))) +^(°> 72^00)))- 

Carrying the same computations as in the geodesic case, we have 


71,72 


0 ) 


—2Re 


V00) - ^72^)00) -rVfl/ 


We will show that <3?^ is non identically vanishing for generic £. 
Evaluating the above formula at r — 0 yields 


t;„ 72 (0) = —2Re (f mg-Wt - vLr'Q) • 
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Since we are assuming 71 7 ^ 72 we can use the exact same ideas as 
before to show that 

for a set of £ with full measure in the discontinuity set. Being a real- 
analytic, non identically vanishing function, < 3>1 ?| 72 (r) has a holomor- 
phic extension to an open complex domain 

[-r 0 , +r 0 ] CllcC 


and by further shrinking hi we can assume that it has finitely many 
zeros Zi,... ,z,j ( repeated with multiplicity ) in hi. The map 

„ ^ 71 , 72 ^) 

nUi*-*) 

is holomorphic, non vanishing on hi and therefore there exists C > 0 
such that for all r G [—To, +fq], 


KM>c 


n<* - z ^> 

3 = 1 


We can then apply Lemma T 6 to show that for generic £, all 71 7 ^ 72 

lim J 7 l) 72 (A)=0. 

A—>-+oo 


We now need to treat the second phase function © 71 ) 72 (V). Performing 
similar calculations we have 


= 2Re 



l^(r)| 2 (£i + £ 2 ) - 2 t{r) - 2 i(r)^ 2 + £1 + £2 
(1 -|£(r)P)(£(r)-£ 1 )(£(r)-£ 2 ) 


where we have set for simplicity £1 = g X 7 X x £, £2 = g x y 2 x £. We 
obtain for r = 0 , 


e;„ 7a (0) = 2Re (nox-rS-y + <r V?)) • 

We want to show once again that for a generic choice of £, this is not 
0. First, remark that we cannot have for all £ G S 1 


9 = -9 X 7i *£• 


Indeed, such an identity would imply (by analytic continuation) that 
for all z G HI 2 , 

W 1 7 2 " 1 7W(-) = ~z. 

If 7 ! = 72 we clearly have a contradiction while if 71 7 ^ 72 this formula 
would show that 7^71 is an elliptic isometry , simply because it is 
conjugated to * G — z, which is elliptic. Because T is a convex co¬ 
compact group whose elements are all hyperbolic (except identity), we 
have again a contradiction. Therefore 

f/ _1 72 _1 £ = -^ _1 7 i _1 £ 
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can hold for at most two points in S 1 . By removing a countable set of 
the discontinuity domain, we can rule out this case. We are left with 
the case 

which can be treated as in the previous section by using an orientation 
argument. Discarding another countable set of points, we can make 
sure that for all 71,72 € T, 


e' (o) / 0. 


We can then use the same arguments as before and apply Lemma 3.6 
to get decay of oscillatory integrals 


lim J 7 l , 72 (A) = 0. 


A—^-l-oo 


To conclude this section we point out that it is unclear to us whether 
Proposition |3.1| holds for general analytic curves, which prevents us 


from extending the lower bound of Theorem 1.1 to analytic curves 


However as pointed out in the next section, it works without major 
modification for the upper bound, extending the upper bound of The¬ 


orem 1.1 to real-analytic curves. 


4. Counting intersections of nodal lines with geodesics 

In this section we prove Theorem |1.1[ using the previous equidis- 
tribution result. We assume that C = g([— 1, +1]) is a fixed geodesic 
satisfying £-NS, and that <5(T) < ^ as we did before. 


4.1. The lower bound. Let Co C C be a geodesic segment given by 
Co = g([— r o, + r o])- We pick J C [—r 0 , +r 0 ] so that the conclusion of 
Proposition 3.1| holds. Since we have 

j \ F \(g(r),Z)\dr > \\Fx\\ll {g{J)) J^F x (g(r),£)) 2 dr, 

remembering the bound (|2| we can use Theorem 3. 5| which says that 

lim f (F x (g(r),£)) 2 dr = \ f E^gir^Qdr, 

A ^°° J j J j 

to conclude that one can find C > 0 such that for all A large, 


\F\(g(r),£)\dr > C. 

Let iV(A) > 0 be the number of zeros of r t —> F\(g(r), £) in the interval 
Int(J). By writing 

A(A) 

1=0 
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where r H- F\(g(r), £) has constant sign on each Ja, we deduce by 
Proposition 3.1| that 

r r 

o <c< I F x (g(r),Z)\dr=Yl / (g(r),£)dr 

Jj i =0 J Jt 


rP 


<(N(X) + 1) sup / F\(g(r),£)dr 

a</3GJ Ja 

which implies that for all A large enough, 

N{ A) > C' A, 

and the proof of the lower bound is done. 


< 


C(N( A) + 1) 
A 


4.2. The upper bound. As we said in the introduction, we will need 
to use analyticity to prove the upper bound on the number of intersec¬ 
tion of nodal lines with geodesics. We will therefore start by proving 
the following fact, which is a way to ’’complexify” restrictions of eigen¬ 
functions F\ to geodesics. 

Proposition 4.1. Let G = g([— 1, +1]) be a geodesic. Then for all A G 
M, the map z H> F\(g(z), £), defined on (—1, +1), admits a holomorphic 
extension to the unit disc D, which is denoted by F\ }g (z,£,). Moreover, 
for all compact subset K C D, there exist fix , Ck > 0 such that for all 
A > 0, we have 

sup |F Aifl (z,OI < C K e pKX . 

z£K 

Proof. We recall that for all r G (—1,+1), we have the convergent 
series expansion 

fA(s(rU) = J]eh' (0 ^ w) co S (AB £ (0,79( r ))). 

7er 

Since we have 


^(0,7 g(r)) = Bfi0,jg(0)) + B g -i 7 -i ( (0,r), 
it is enough to continue analytically 

r kg = log 


V ~ 0 _1 7 _1 f| 2 , 

We set for simplicity rj := g l 7 _1 £ and for all z G (—1, +1), 

1 — zfi 1 — z 2 

\z — 7]\ 2 (z — rj)(z — rj )' 

Clearly G v (z) extends holomorphically to the unit disc D, where it does 
not vanish. We can therefore define a complex logarithm by setting for 
all :GD 


( 8 ) 


~h(G v )(z) := 


G' V {Q 

G V (C) 


d( = z 


ml G' v (zt) 
) G v (zt) 


dt. 
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We obtain a holomorphic function L (G v )(z) on D which has the fol¬ 
lowing properties: 

• V r £ (-1, +1), L (G v )(r) = log G v (r) = B n { 0,r). 

• Vz £ O, e L ( Gr, )( z ) = G v (z). 

By using formula (J 8 |, one can check that for all 0 < 77 < 1, 

sup |L(G r? )(^)| < C(r i), 

\z\<ri 

where C(r i) is uniform in r) := g l 7 -1 £. Writing 
cos (A^( 0 , 75 f( 0 )) + AL(G, ? )(^)) 


= cos (AOg(0, 7 ^( 0 ))) cos(AL(G' r? )(z))-sin(AS 5 ( 0 , 75 f( 0 ))) sin(AL(G , r? )(^)), 
and using the bounds for all z £ C, 

| cos(^)| < 2 e |Im(z)l , | sin(^)| < 2 e |Im(z)l , 
we deduce that for all \z\ < ry and A > 0, 


cos (AO ? (0,70(0)) + AL(G„)(z)) I < C(n)e^\ 


Combining this last bound with Lemma 2. 1| shows uniform convergence 

{\z\ < ri} 


on 


of the above series, hence holomorphy and the claimed bound. □ 
Notice that for a more general real-analytic curve, a similar state¬ 
ment follows straightforwardly, with the difference that it will hold on 
a smaller domain 0 C D. 


The rest of the proof of the upper bound on the number of intersections 
of nodal lines with Co C C will follow from Theorem 3J3 combined with 
Jensen’s formula. The version of Jensen’s formula we will use is the 
following. 


Proposition 4.2. Let f be a holomorphic function on the open disc 
D(w,R), and assume that f(w ) ^ 0. let Nf(r) denote the number of 
zeros of f in the closed disc D(w,r). For all r < r < R, we have 

Nfi ?) < \ og l {r /r) (JfrJ 0 l °g\f(™ + re ie )\dO-\og\f(w)\^ . 

For a reference on Jensen’s formula, we refer the reader to the clas¬ 
sics, for example Titchmarsh [IB]. Let C 0 = g([—r 0 , +r 0 ]) be a geodesic 
segment as above, with 0 < ro < 1. Fix e > 0 so small that r 0 + 3e < 1 
and set 

ri = r 0 + e, r 2 = r 0 + 2e, r 3 = r 0 + 3e. 

If D(w,r ) denotes the complex open disc with center w and radius r, 
we then have for all x £ [—e, +e] 


.D(0, r 0 ) C D(x,r i) C D(x, r 2 ) C 0(0, r 3 ) C D. 
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Let N( A) denote the number of zeros of r i—>• F\(g(r),£) in the interval 
[—r 0 , +r 0 ]. By applying Theorem 3.5 on the short interval [—e, +e], we 
have 




f+e 


lim 

A—^+oo 


(F\(g(r),€)) 2 dr = \ j E x {g{r),i)dr, 
which shows that for all A large enough we have 

0 < C e-=l(^~ [ E 1 (g(r),^)dr\ < sup \F x (g(r), £)\. 

\ Ze J-e J re[-e,+e] 

For all A large, we denote by X\ E [—e, +e] a point such that 

\F\(g(x x ),0\ = sup |F A (c/(r),OI- 

re[— 6,+e] 

Applying Jensen’s formula to F\ g (z,£,) on D(x\,ri) C D(x\,r 2 ), we 
have 

N(x) 5 MG7G) (-hi io s\h,^ + r 2 e“,md») 

(log \F\(g(x\), £) |) 


log(r 2 /ri) 


< 


log(r 2 /ri) \ M< r3 


sup log \F X: g(z, ^)| +log(C, 




Using the estimate of Proposition |4.1[ we then deduce that as A goes 
to infinity, 

N{\) = 0{\), 

and the proof is completed. □ 


To deal with more general real-analytic curves which extend holo- 
morphically to a smaller domain C D, we just need to replicate the 
same argument with several discs instead of a single one. We omit it 
for simplicity. 


5. Counting nodal domains 

5.1. The non-elementary case. Let us introduce some notations. 
We assume in this section that T is non-elementary. We will work on 
the universal cover HI 2 , so that the convex core A" 0 is the image under 
the covering HI 2 —> r\HI 2 of a compact geodesic polygon 

7 C H 2 . 

The polygon 7 has finitely many sides which are geodesic segments, 
see the picture below for an example such a polygon in HI 2 = D, the 
gray hyperbolic octogon is 7. 
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We choose £ G S 1 \ A(r) such that the upper bound of Theorem 


1.1 


is valid on the full boundary 07, which can be done for a set of full 
measure. We recall that the nodal domains of F\(z, £) : HI 2 —>■ M are 
by definition the connected components of 

M 2 \{F x (z,0 = 0}. 

The nodal domains 7) which do intersect 7 fall into two categories. 
Either 

V n 07 ± 0, 

and thanks to Theorem |1.1| there are at most 0(A) of them, or we have 

V C Int(fP). 

In that case, since F\ has constant sign on 7), the eigenvalue 

/i = 1/4 + A 2 


must be the first eigenvalue of the hyperbolic Laplacian Ajp on 7) for 
the Dirichlet boundary problem. 


f A H 2^ = 

['0 = 0 on 07). 

Let Ai(D) denote the smallest eigenvalue for the above Dirichlet prob¬ 
lem. We will use the following key lower bound. 


Proposition 5.1. Fix eo > 0, then there exists Cq > 0 such that for 
all domain SlcB 2 with Vol(D) < eo 


Ai(fi) > 


Co 

Vol(fi)' 
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Proof. We first use the Faber-Krahn inequality for domains in H 2 , see 
Chavel jjj p. 87. It is valid on simply connected spaces of constant 
curvature. If 17 is a compact domain of H 2 with piecewise C°° bound¬ 
ary, then the first Dirichlet eigenvalue Ai(f2) of the Laplacian satisfies 
Ai(fl) > Ai(-D), where D is a geodesic disc with same (hyperbolic) vol¬ 
ume. The game is now to prove a lower bound for the first eigenvalue 
on a disc D of the hyperbolic plane, for small values of the radius. We 
use the disc model for HI 2 and can assume that D = D(0,r) (euclidean 
disc) is centered at 0. By the min-max principle, we have 


Ai(-D(0, r)) 


• nf f D ^(A H 2(^)dVol 

<^o eCg°(D) J D tp 2 dVo\ 


But we have 


( p(A^2(p)dVol = / (p(A(p)dm , 


' d 


id 


where m is the Lebesgue measure and A the positive euclidean Lapla¬ 
cian, while 


Lp 2 dV ol = / ip 2 {z) 


4 dm(z) 


< 


JD JD (l - \z\ 2 ) 2 ~ (1 - rl) 2 

as long as r < ro < 1. We therefore have 


tp 2 (z)dm(z), 


i D 


Ai(O(0,r))> (1 4 r ° )2 Ar(g(0.r)), 

where A® uc denotes the first Dirichlet eigenvalue for the euclidean Lapla¬ 
cian. A simple change of coordinates in the min-max then shows that 


\? c (D(0,r))> 


Ar(D(0,l)), 


Using the formula for the hyperbolic area of D(0, r) 
Vol(fi) = Vol(D(0 ,r)) = l;/ ” 


1 — r 2 


shows that 


Ai(fi) > 


Ml - rg)Ur(P(0,1)) 

Vol(fi) 


and the claim is proved. □ 


Going back to the proof of the upper bound, let (Dj) ie / be the (fi¬ 
nite) collection of nodal domains T>i that are inside Int(CP). By volume 
comparison, we have 

Vol(U ie /Di) = Y, Vol(Di) < Vol(CP). 
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Let J C / be the set of indexes such that for all j G J, Vol(Dj) < eo- 
By Proposition 5.1 we get 


#(J)Co 


< Vol(fP), 


1/4 +A 2 

which obviously shows that #(J) = 0(A 2 ). Similarly we have 
#(/\J)<eo 1 Vol(3>) = 0(l). 


As a conclusion we have shown that the total number of nodal domains 
that intersect 7 is 0(A 2 ), thus completing the proof of the upper bound. 


5.2. The cylinder case. Here we assume that T is an elementary 
group so that X = r\HI 2 is a hyperbolic cylinder. We denote by C 0 the 
unique closed geodesic on A". Fix r > 0. The collar 

C(r) := {z G X : dist(+ C 0 ) < r} 

is the image under the projection n : H 2 G X of a domain 7 in H 2 
whose boundary is piecewise circular (not totally geodesic). 



More precisely, we can (up to a conjugation by an isometry) assume 
that Co lifts in HI 2 to the segment (—1,+1), so that by a classical 
formula (see Beardon [lj, p.163) we have 

dist(+ C 0 ) < r <G> -—r/ < sinh(r). 

1 — \z\ 2 

If Co is the axis of a hyperbolic isometry y 0 and V is the group generated 
by 70 , then a fundamental domain for the action of T is provided by the 
domain of HI 2 which is outside the isometric circles of 70 and 7 0 ~ 1 . The 
domain 7 is then the grey region depicted in the previous picture, which 
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correspond to the intersection of the collar (in H 2 ) with a fundamental 
domain. Since d {P is piecewise real analytic, we can adopt the exact 
same strategy as in the previous proof, by choosing £ such that Theorem 
|1.1| applies, and by arguing the same way, depending on the type of 
nodal domain. 


5.3. Lower bounds and open questions. The first remark that we 
have in mind is that by adapting straightforwardly the combinatorial 
arguments used in Eaim ilO j we can obtain a lower bound for the 
number of connected components of A"o \ for generic £, which says 
that for large A, 


M s (A) > C~ l A. 


Clearly the main input here is the lower bound given by Theorem |1.1| 
and the graph theoretic arguments from in. which are a generaliza¬ 
tion of the more elementary ideas pioneered in [7j. However that kind 
of lower bound is rather irrelevant, because we cannot rule out the 
fact that these connected components could very well come from a sin¬ 
gle nodal domain which would intersect several times the convex core 
Xq. These issues are already present on compact manifolds where one 
has either to use symmetries or boundary conditions to rule out these 
pathologies. 


From the numerics one can formulate the following list of open ques¬ 
tions which seem to be relevant. 

• It seems that for compact sets K with non empty interior which 
are in the vicinity of £, the number Mk{ A) of nodal domains 
that intersect K obeys the growth rate Mk( A) x A. 

• Is the number of compact nodal domains finite ? Do compact 
nodal sets remain in a compact part of the surface, uniformly 
in A ? 

• Prove that there exist compact nodal domains, if A is large 
enough, start with the elementary group case. 

• Prove or disprove that the number of compact nodal domains 
inside the convex core is, as A —* +oo, greater than CA 2 , for 
some C > 0. This question could be tested numerically. 

• In the plot below we have found for A = 150 some examples of 
non-simply connected compact nodal domains. Can the topol¬ 
ogy be arbitrary ? 
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